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ABSTRACT 


We  develop  a  new  method  for  the  calculation  of  the  atom 
scatterins  analog  of  the  Debye-Waller  factor.  Unlike  X-ray  and 
neutron  scattering  the  properties  of  an  atom  scattered  by  a  solid 
surface  cannot  be  computed  by  perturbation  theory,  therefore  the 
simple  Debye-Wal ler  theory  cannot  be  applied.  Nevertheless  the 
Oebye-Waller  phenomenon  has  a  close  analog:  the  elastic 
intensity  is  depressed  due  to  the  uncorrelated  part  of  the 
thermal  motion  of  the  lattice  atoms.  To  compute  this  effect  we 
develop  a  time  dependent  scattering  theory  in  which  the  quantum 
properties  of  the  scattered  atom  are  described  by  propagating 
coherently  an  ensemble  of  wave  packets  and  lattice  motion  is 
simulated  by  a  classical  Langevin  equation.  Applications  are 
made  to  He  and  Ne  scattering  from  a  surface  whose  lattice 
dynamics  mimics  that  of  Pt(lll)  but  whose  corrugation  was 
slightly  increased  to  enrich  the  diffraction  structure. 


A 


Introduction 


I . 

The  diffraction  of  atone  by  solids^  has  many  useful 

applications  in  surface  science,  particularly  as  a  probe  of  both 

surface  structure  and  dynamics.  The  shapes  and  intensities  of 

the  diffraction  peeks  are  very  sensitive  to  the  tbernal  notion  of 

the  surface  atons.  In  this  paper  we  present  a  theory  for  the 

effect  of  this  notion  on  the  elastic  diffraction  peak 

probabilities;  that  is,  for  the  aton  scatterine  analog  of  the 

Debye-Waller  factor  appearlne  in  x-ray  and  neutron  scattering. 

2 

Prior  work  has  discussed  such  effects  by  taking  over,  with 

alight  nodif ications  the  Debye-Waller  fornula  fron  neutron 

2f 

scattering  theory.  Weinberg  ,  Hoinkes,  Nahr,  and 
3  2 

Wllsch  ,  and  others  used  the  resulting  equation  to  extract 
values  for  the  surface  Debye  tenperature  of  the  lattice  atons  and 
the  results  obtained  are  reasonable.  So,  while  this  procedure  is 
not  terribly  nisleading  it  is  suspect,  fron  a  theoretical  point 
of  view,  on  several  counts.  The  derivation  of  the  Debye-Waller 

4 

factor  in  neutron  scattering  relies  ..eavily  on  the  fact  that  the 
neutron-nucleus  interaction  can  be  described  by  the  Feral 
pseudopotential .  ^  As  a  result  the  collision  of  the  neutron  with 
a  given  nucleus  can  be  calculated  by  using  the  Golden  Rule  and 
the  potential  has  zero  range.  Furthernore,  the  neutron" lattice 
interaction  is  so  weak  that  we  can  assume  that  the  neutrons 
passing  through  the  sample  have  been  scattered  once  (i.e.  by  one 
nucleus)  or  not  at  all.  Clearly  none  of  these  conditions  hold 
for  atom  scattering.  A  nore  detailed  analysis  of  the  difference 
between  aton  and  neutron  scattering  and  of  the  reasons  why  such 
differences  affect  the  Debye-Waller  factor  can  be  found  in 
Appendix  1 . 

a 

The  application  of  the  usual  coupled  channels  method  to 
collisions  with  a  surface  undergoing  thernal  notion  seems  im¬ 
practical  unless  the  phonon-particle  interactions  are  treated 
within  perturbation  theory  (the  distorted  wave  approxina- 


tion).  It  is  not  clear,  however,  whether  perturbation  theory 
is  adequate  for  this  situation.  One  would  therefore  like  to  have 
a  non-per turbative  approach  to  this  problem.  One  possibility  is 
to  treat  the  motion  of  the  lattice  classically  and  the  scatterer 
by  quantum  theory.  In  general  such  an  approach  has  difficulties 
since  a  replacement  of  some  of  the  quantum  variables  by  their 
classical  trajectories  leads  to-abieuities  in  defining  the  proper 

O 

classical  (Newton)  equations  :  the  theory  does  not  (in  the 
simplest  cases)  provide  a  feedback  mechanism  by  which  the 
classical  degrees  of  freedom  are  ’’informed"  that  they  have  to 
change  their  energy  when  the  quantum  ones  are  changing  theirs. 

One  way  of  avoiding  such  difficulties  is  to  treat  the 
quantum  degrees  of  freedom  by  a  method  which  is  very  compatible 
with  classical  mechanics.  The  possibility  which  we  pursue  here 

Q 

is  the  wave  packet  method  developed  by  Heller  and  his  coworkers 

and  applied  to  diffraction  by  a  static  lattice  by  Droshlagen  and 

Heller  (DH).^®  The  method  consists  of  writing  the  Incident 

planar  wave  as  a  sum  of  Gaussian  w.ave  packets,  propagating  them 

independently  by  using  the  time  dependent  Schroedinger  equation 

and  constructing  the  scattered  wave  function  by  adding  coherently 

the  scattered  packets.  The  theory  reduces  the  wave  packet 

propagation  to  the  calculation  of  the  time  evolution  of  the 

position  of  the  center  of  the  wave  packet,  its  mean  momentum,  and 

its  width  and  phase.  The  center  of  the  packet  and  the  momentum 

evolve  according  to  the  classical  equations  of  motion,  the  phase 

at  the  end  of  the  trajectory  is  the  classical  action  and  the 

width  satisfies  a  first  order  differential  equation.  All  these 

quantities  depend  on  the  instantaneous  position  of  the  lattice 

atoms  through  the  classical  atom-lattice  interaction.  We 

propagate  the-  lattice  atom  positions  by  using  the  Adelman-Dol 1 - 
1 2 

Tully  method  (ADT)  which  provides  a  Langevin  equation  for  the 
motion  of  the  surface  atoms  involved  in  the  collision  with  the 
probe  atom.  As  a  result  of  combining  the  Langevin  equation  for 
the  surface  atoms  with  the  wave  packet  description  of  the 


Incident  particle  the  motion  of  the  packet  is  modified  as 
follows:  the  position  and  the  momentum  p^  of  its  center  are 

now  stochastic  variables  evolving  according  to  a  classical 
Langevin  equation;  the  phase  which  depends  on  r^  and  p^  is  also 
randomized,  thus  creating  a  dephaslng  which  leads  to  a  decrease 
in  the  diffraction  peaks  Intensity  and  an  increase  in  the 
background  scattering.  The  wldth^which  provides  a  space 
dependent  phase  and  an  amplitude-  also  becomes  a  random  variable 


THE  METHOD 


H  . 


Since  we  plan  to  treat  the  lattice  motion  classically  we 
start  with  the  customary  classical  trajectory  approximation 
method  (CTA).  If  the  Hamiltonian  depends  on  two  sets  of 
variables  x  and  y  and  we  want  to  treat  y  classically  then  we  use 
the  "reduced"  Hamiltonian 

H(x.y(t))  -  -{h^/2m)  7^  ^  V(x,y{t))  .  (II. 1) 

2  2 

The  kinetic  energy  operator  -(h  /2m)7y  is  removed  from  the 
Hamiltonian  since  this  quantity  is  included  in  the  classical 
equation  of  motion  for  y(t).  The  potential  energy  and  the  wave 
function  i)i(x;y(t))  depends  on  the  classical  variable  y(t) 
parametrically.  In  the  present  case  the  quantum  variable  x  is 
the  atom  position  ?  and  the  group  of  classical  variables  y(t)  is 
the  set  R ^ ( t )  of  lattice  atom  positions. 

Following  DH  we  write  the  time  dependent  wave  function 
i|)(7,t)  of  the  incident  atom  in  the  form 

»li(T;t)  -  Z  G^(r:  {X^(t)))  (II. 2) 

tt  0( 

a 

where  G  is  a  Gaussian  wave  packet  (GWP)  given  by 
oc 

G  (T.iX  (t)})  -  exp{  { 1/h)  [  (T-lr  ( t )  )  A  *  ( t )  •  (^-T  (  t )  )  - 

a  ft  ft  ft  ft 

p  ( t ) • (7-7  ( t ) )  +  r„(t)]}  (1 1. 3) 

ft  ft  ft 

The  set  (X  (t))  is  a  symbolic  notation  for  the  time  dependent 

parameters  A  ,  r  (t),  p  (t)  and  y  (t).  The  parameters  r  (t)  and 
^  0(  a  oe  oc  a 

p  (t)  are  the  expectation  values  of  the  position  and  momentum  of 

ft 

the  packet 


-  fdT  G^(‘?:{X  (t)})* 

V  ft  ft 

p„(t) 

) 

; 

G„(r:{X  (t)} 
a  oe 


(  II  .4) 


The  complex  quantity  y  (t)  generates  the  phase  Rey  ( t ) /h  and  the 

ft  ft 

amplitude  exp{-Iffl  y  (t)/}i}  of  the  packet.  The  three  dimensional 
square  matrix  A^(t)  controls  the  width  of  the  packet  (through  its 
imaginary  part),  generates  a  space  dependent  phase  (through  its 
real  part)  and  couples  the  variables  x.y,  z  (in  the  interaction 
region)  through  its  off-diagonal  elements. 


The  time  evolution  of  the  wave  function  ii>(r;t)  is  obtained 
by  what  we  call  the  simplest  Keller  method  (SHM).  This  forces 
the  parameters  to  vary  in  time  so  that  each  Gaussian 

satisfies  the  time  dependent  Schrodinger  equation 

[-(h^/2m)  7^  +  V(r:{^^(t))]  6^( r ; {X^( t ) ) )  -  ih  G^(r:{X^})  . 

r 

(  II .5) 

We  use  the  CTA  Hamiltonian  in  which  the  lattice  atoms  positions 
^^(t)  are  treated  as  classical  fields. 


To  obtain  equations  for  the  tine  evolution  of  X^(t)  SHM 
expands  V(?;(R^(t)})  in  a  power  series  around  the  instantaneous 
packet  position  r  (t)  and  neglects  the  terns  higher  than  second 

ft 

order  in  (r  -  r  (t)).  By  introducing  this  expansion  in  (II. 5) 

QC 

and  equating  the  coefficients  of  (r  -  r  (t))",  n  =  0,1,2  to  zero 

oc 

SHM  gives 


r^(t) 


p^( t ) /m 


(11.5a) 


-3V(r  (t):{R.(t)})/3r^(t) 

QC  I  OC 


( II .5b) 


P„{t)  - 


A  (t) 
a 


and 


-(2/»)  A^(t)-  A^(t) 

a  OC 

-  (1/2)— 2 - -  .V(T  {t):{R.(t)}) 

3?^{t)  3T  (t)  “  ^ 

ft  ft 


(11.5c) 


y^(t)  -  (lh/m)Tr{'^(t)}  +  (l/2m)t  (t)^  -  V  (7  ( t )  :  (R,  ( t ) }  )  . 

CC  OC  OC  OC  1 


( II  . 5d) 

2 

Here  3  V( r_; (R . } )/3r  3r  Is  a  matrix  whose  xy  element  is 

2  ft  1  ft  ft 

3  V/ 3x  ( t ) 3y  ( t ) ,  where  x  (t).  y^(t)  and  z  (t)  are  the  components 

ft  ft  ^  ft  ft  ft 

of  the  vector  r  (t). 

ft 


These  equations  indicate  that  r  (t)  and  p  (t)  follow 

ft  ft 

classical  trajectories  with  the  potential  V(r  (t);  (R.(t)}). 

OC  1 

Furthermore,  one  can  show  that  the  phase  y  (t)  is  equal  to  the 
classical  action  along  these  trajectories.  The  width  “T  (t)  has 

ft 

its  own  dynamics,  with  no  classical  analog,  and  depends  on  the 
force  constant  tensor  (i.e.  the  second  order  derivatives  of  the 
classical  potential). 


To  solve  the  equations  (11.5)  we  must  provide  a 

prescription  for  the  classical  motion  of  the  surface  atoms.  For 

this  we  use  a  method  proposed  by  Adelman  and  Doll^^  and  developed 
1 2 

by  Tully.  This  provides  a  Langevin  equation  for  the  motion  of 
the  lattice  atoms  suffering  the  brunt  of  the  collision. 


III.  LIMITATIOHS 


In  order  to  see  Laue  diffraction,  our  state  must  be 
sufficiently  coherent  in  x  and  y  over  several  unit  cells  to  probe 
the  surface  periodicity.  We  accomplish  this  by  overlapping  many 
packets  in  a  two  dimensional  grid  covering  a  single  unit  cell, 
choosing  their  initial  phases  such  that  together  they  simulate  a 
plane  wave  and  u-sing  Bloch's  theorem  to  extend  this  coherence 
over  many  unit  cells.  Following  DH ,  this  is  done  by  making  use 
of  the  identity: 

1^0 

4)  (  r  ,  0  )  =  ce 

«  c'Jd?  exp{  l/h[  dr-r”  ) '  ( r -?  )  +  ik  •  {T-T  ) 

0  O  O  O  0  o 

*  (III.l) 

0  o 

Replacing  this  integral  by  a  sum  over  many  equally  spaced  points 
r^  produces  a  grid  of  packets  which  mimic  a  plane  wave.  These 
packets  are  scattered  from  the  surface  and  the  final  state  ij)  ( r , 
t)  is  projected  onto  various  outgoing  plane  wave  states  X(r), 
where 

X(r)  =  e^^'^  (  III  .2) 


Thus,  the  probability  of  being  in  a  plane  wave  state  having  some 
wave  vector  is  the  square  of  the  S  matrix 


P(q)  •  |S(q)l2 

where 


(III  .3) 


S(q)  -  <X|4l> 


(III  .4) 


One  can  show  that  for  the  case  of  a  static  lattice,  it  is 
sufficient  to  scatter  from  a  single  unit  cell  in  order  to 
calculate  the  diffraction  probabilities.^®  In  this  paper,  we 
also  use  a  single  unit  cell,  that  is,  we  use  the  ADT  method  to 
propagate  a  four  atom  primary  zone  on  the  surface.  The  result  is 
that  we  suppress  the  long  range  correlations  between  the  motions 
of  the  surface  atoms,  and  we  therefore  cannot  reproduce  the 
effect  of  thermal  motion  on  the  diffraction  line  widths. 

However,  this  does  not  affect  the  "Debye-Waller"  factor  for  the 
elastic  diffraction  peaks,  since  this  arises  from  the 
uncorrelated  motion  of  individual  atoms.  More  is  said  about  this 
in  Appendix  2. 

Two  other  limitations,  should  be  mentioned:  (1)  the  use  of 

a  single  layer  of  packets  along  the  z  coordinate  makes  the 

incident  state  partially  coherent.  The  manner  in  which  this 

1 3 

affects  diffraction  was  discussed  in  a  previous  paper.  (2)  the 

13  14 

use  of  SHM  nay  introduce  some  errors;  *  however*  existing 

10a 

calculations  agree  with  some  exact  results. 


10 


IV.  numerical  results 

The  procedure  of  computation  is  as  follows: 

1.  Makine  use  of  Eq .  III.I,  we  write  our  initial  plane 
wave  state  as  a  grid  of  packets  covering  the  area  of  one  unit 
cell.  A  15  X  15  array  is  usually  sufficient.  The  distance 
between  the  center  of  a  packet  and  the  surface  is  6  A  or  more. 

2.  Each  packet  is  then  propagated  individually,  is 
scattered  from  th^  primary  zone,  and  its  contribution  to  the  S 
matrix  (Eq.  III. 4)  is  computed.  The  initial  conditions  and  the 
random  fluctuations  which  determine  the  notion  of  the  primary 
zone  atoms  are  the  same  for  each  packet. 

3.  The  S  matrix  contributions  for  all  trajectories  are 
added  (coherently)  and  the  result  is  squared  to  find  the 
scattering  probabilities. 

4.  The  procedure  is  then  repeated. several  times,  using 
each  time  a  different  set  of  initial  conditions  and  fluctuations 
for  the  lattice  atoms.  The  results  are  averaged  until 
convergence  is  seen.  This  requires  only  a  few  full  runs  below 
100  K,  and  as  many  as  30  or  more  above  300  K. 

Using  this  procedure  we  examine  the  diffraction  of  He  and 
Ne  from  a  Pt(lll)  surface.  The  atom-surface  potential  consists 
of  two  terms ; 


V(r) 


S  [ 
i 


r-R. 


12 


(IV.l) 


The  first  is  a  sum  of  individual  Leonard- Jones  6/12  interactions 
with  several  surface  atoms  located  at  R^.  Four  of  these  are  the 


11 


primary  zone  atoms,  and  the  rest  (usually  10  are  sufficient) 
surround  the  primary  zone  to  help  define  the  surface  corrugation. 

12ft 

As  discussed  in  Tully's  paper,  the  primary  zone  is  moved  along 

the  surface  so  as  to  be  always  below  the  Impinging  packet.  The 

second  part  of  Equation  (IV. 1)  Is  a  bulk  interaction,  resulting 

from  a  pairwise  summation  of  Lennard-Jones  6/12  interactions  with 

all  the  atoms  in  the  solid.  It  is  only  a  function  of  z,  and  has 

the  proper  z  long  range  attractive  part.  The  parameters  c^ 

through  c.  were  chosen  to  produce  a  well  depth  of  .125  kcal/aole 
4 

and  a  corrugation  of  .19  A  for  a  Ne  bean  energy  of  1.0  kcal/mole. 

This  particular  fora  for  the  potential  was  chosen  not  to  mimic 

the  actual  He/Pt(lll)  system,  which  has  negligible  off  specular 
1 5 

scattering,  but  to  provide  a  model  system  with  reasonable 

corrugation  for  study.  The  primary  zone  uses  the  parameters 

12 

calculated  by  Tully  to  simulate  the  actual  Pt(lll)  surface. 
Unless  otherwise  noted,  the  beam  energy  was  taken  to  be  1.0 
kcal/aole,  with  an  incident  polar  angle  e  -  30°  along  the 
negative  x  direction  (^  -  0). 

When  we  calculate  S(q)  we  generally  choose  q  and  q  at  or 

*  y 

near  a  diffraction  peak,  fixing  q  by  energy  conservation.  We 

z 

could  similarly  choose  a  detector  angle  (e,^)  and  fix  the 
magnitude  of  q,  that  is  |q|  -  I  '  **kere  |  is  the  magnitude 
of  the  initial  wave  vector.  This  is  only  valid  in  the  limit  of 
elastic  scattering.  If  we  want  to  describe  experiments  in  which 
the  energy  of  the  scattered  beam  is  not  analysed  we  must 
integrate  over  all  final  |q|'s,  allowing  for  all  possible  final 
He  energies.  Thus  the  detection  probability  at  some  angle  (e,^) 

1  s 


P(«,^)  -  :  d|q|  q^  P(|q|.«,«)  (IV. 2) 

0 

The  easiest  way  to  determine  this  would  be  to  compute  P  at 
some  (e,^)  for  several  jqt's,  and  add  then  up  with  the  proper 
weighting.  We  again  encounter  a  problem  with  beam  coherence. 


Our  initial  state  has  a  finite  width  in  z  and  covers  only  a 
single  unit  cell  in  x  and  y.  Due  to  this  finite  spatial  extent, 
its  nonentuB  is  not  so  sharply  defined  as  in  a  planar  wave.  Thus 
the  final  distribution  of  probabilities  as  a  function  of  |'q| 
depends,  upon  how  Buch  area  our  initial  state  covers.  Consider 
Figure  1  where  we  plot  P(|q|.o.^)  for  various  values  of  |q| 
centered  about  |q|  -  |  and  (e.^)  corresponding  to  the  specular 

peak.  As  we  Increase  N.  where  i|>(r.o)  now  covers  NxN  unit  cells, 
the  distribution  sharpens  up  as  expected.  Thus  P(d,^)  in 
equation  (IV. 2)  is  a  function  o(  N.  Since  it  Bight  require  very 
large  N's  to  Bake  this  dependence  go  away,  we  choose  |*q  |  >  l^g  I  * 
What  we  calculate  then  is  the  elastic  DW  effect,  which  is 
equivalent  to  doing  a  tlae  of  flight  diffraction  experlBent.  We 
conflrB  that  it  is  N  that  doBinates  P(iq|)  (for  SBall  N)  in 
Figure  2.  We  plot  two  distributions  siailar  to  figure  1  for  the 
2x2  cell  case,  but  at  different  tesperatures .  Norsally  the  400 
K  case  would  have  a  broader  distribution,  but  both  have  the  sane 
N  -  2  Halting  width. 

In  teras  of  the  diffraction  peak  widths  in  real  or  k-space, 
we  see  the  saae  behavior.  In  Figure  3  we  plot  P  vs  AK  in  the 
vicinity  of  the  (1,0)  peak.  Alt  is  the  change  in  the  parallel 
coaponent  of  the  beaa  aoaentuB,  in  units  of  t,  (the  reciprocal 
lattice  vector  in  the  (1,0)  direction).  We  plot  P  for  both  N-2 
and  3  and  T-lOK  and  200  K.  For  each  N  we  see  the  saae  peak 
intensities  and  DW  effects,  with  only  the  width  being  different. 
Again  the  bean  coherence  doalnates  and  we  find  their  widths  to  be 
t/N,  as  in  Equation  A. 4. 

For  the  above  He  systen,  the  (1,0),  (0,1)  and  (1,1)  peaks 

all  have  a  reasonable  probability  and  we  consider  their  intensity 

vs.  teaperature  in  Fig.  4  where  P(T)/P(0)  is  the  ratio  of  the 

peak  probability  at  soae  teaperature  T,  to  that  at  T-0.  The 

*  T 

curves  reseable  the  fora  e  as  predicted  by  first  order 
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2 

perturbation  theory..  This  ib  not  surprising  since  He  is  a  very 

light  scatterer  and  the  He/Pt  mass  ratio  is  only  about  .02.  At 

higher  energies  and  for  larger  atomic  masses,  we  should  see  more 

deviation.  In  Fig.  5  we  plot  In  P(T)/P(0)  for  the  (1,0)  and 

(1,1)  peaks.  The  data  resembles  the  straight  line  predicted  by 

2 

the  Debye-Waller  theories  based  on  neutron  scattering.  Note 
that  at  higher  temperatures  we  begin  to  see  deviations  from  a 
linear  behavior.  This  is  exactly  the  behavior  observed 
experimentally,  for  such  systems  as  He/Cu(100).  A  point  to  be 
made  here  concerns  the  two  different  slopes  observed  in  Fig  5  for 
the  (1,0)  and  (1,1)  peak  probabilities.  Since  these  off  specular 
peaks  couple  differently  to  the  x  and  y  components  of  the 
displacement  lT(t)  of  the  surface  atom  (they  have  a  different  Ak^ 
and  Aky) .  they  show  a  different  Debye-Waller  effect.  To  properly 
predict  this  using  the  simple  models  would  require  a  knowledge  of 
the  displacements,  which  are  not  equivalent  for  a  non-cublc 
structure  such  as  Pt(lll). 

In  Fig.  6  and  7  we  plot  similar  data  for  the  case  of  Ne  on 

Pt .  The  beam  energy  is  S.O  kcal/mole,  which  gives  the  Incident 

Ne  atoms  the  same  velocity  as  the  1.0  kcal/mol  He  energy 

considered  previously.  Because  of  the  larger  mass  ratio  (0.1)  we 

-see  much  more  energy  exchange  with  the  lattice.  At  200  X  the 

average  energy  change  is  -.32  kcal/mole  and  AE  >  .72 

ms 

kcal/mole.  For  the  above  He  system,  these  values  were  only  .005 

and  .11  kcal/mole  respectively  at  the  sane  temperature.  In 

addition  to  the  specular  we  Include  data  for  the  (0.2)  and 

(-2,3)  peaks,  all  of  which  show  about  the  same  Intensity.  Again 

we  see  a  different  dependence  on  temperature  for  each  peak.  We 

also  see  a  nearly  linear  behavior  in  Fig.  7  up  to  about  500  K 

2 

when  the  intensity  becomes  larger  than  a  simple  theory  would 
predict.  We  also  see  some  deviations,  for  both  He  and  Ne,  from  a 
linear  behavior  at  lower  temperatures,  but  his  deviation  is 
small . 


The  method  proposed  here  seems  to  be  able  to  account  for 
the  temperature  dependence  of  difraction  peak  intensities  and 
predicts  behavior  similar  to  that  observed  experimentally.  In 
order  to  achieve  a  complete  description  of  diffraction  the  wave 
packet  method  should  be  extended  to  Incorporate  two  effects;  th 
broadenine  of  the  peak  by  energy  and  parallel  momentum  transfer 
and  the  resonances  caused  by  temporary  trapping  at  the  surface 
(selective  adsorption).  As  we  have  already  discussed,  the 
inclusion  of  the  first  effect  is  automatically  achieved  by 
increasing  the  size  of  the  primary  zone  and  by  using  enough 
Gaussians  to  cover  the  whole  zone.  This  leads  to  a  very 
substantial  Increase  in  the  computer  time.  The  Incorporation  of 
selective  adsorption  requires  new  theoretical  developments. 
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APPENDIX  1 


We  discuss  here  the  differences  between  the  present 
approach  and  earlier  work  which  modified  equations  valid  for  X- 
ray  and  neutron  scattering  to  take  into  account  in  an  approximate 
way  certain  features  specific  to  atom  scattering. 


The  simplest  way  to  obtain  the  DW  factor  is  to  examine  the 
scattering  of  a  projectile  by  the  thermal  average  of  the  lattice 
projectile  interaction.  If  we  use  a  pair-wise  potential  V("r)  -  2 
V(7  -  R  )  we  can  write  its  average  value  as 


<V( r  )> 


dk  -  Ik*  r 
- 3  * 


i1c.‘^° 
VCiTxS  e  “ 
a 


e-wdT) 


(A.l) 


We  have  written  here  V(r)  in  terms  of  its  Fourier  transform  V(k) 

and  performed  the  average  by  using  Bloch's  formula  <expQ>  » 

2 

exp[<Q  >/2],  where  Q  is  any  linear  combination  of  the  phonon 
creation  and  anlhllatlon  >.perators.  The  DW  exponent  is 


W(1<)  -  (1/2)  <(1J*6R)^> 

where  6?  is  an  operator  representing  the  displacement  of  the 
lattice  atom  with  respect  to  its  equilibrium  position  (l.e. 
lattice  site). 


For  X-ray  and  neutron  scattering  by  the  bulk  of  a  lattice 
we  can  use  perturbation  theory;  the  transition  amplitude  from  an 
initial  momentum  k^  to  a  final  momentum  is  then  <'k^  |  <V>  |^^  > . 
and  this  leads  to  a  cross  section  proportional  to 

!V(q)i^  exp{-2W(q)  )  IZ  exp{iq*‘^  >|^  6  [  ( /  2m  )  (  k  J-k^  )  ]  ,  (A. 2) 

a 

where  q  s  k^  -  k^  is  the  momentum  transfer:  the  sum  leads  to  the 


Bragg  conditions;  the  5  function  Imposes  energy  conservation. 

One  can  show  that  the  ON  factor  exp{-2W(^)}  will  also 
m,ultiply  the  inelastic  cross  section. 

Physically  the  DW  term  appears  because  the  atoms  are 

displaced  with  respect  to  the  position  of  the  lattice  points 

•*■2. 

according  to  a  Gaussian  distribution  of  width  <6R  >;  the 
displacements  of  any  two  atoms  are  uncorrelated  and  for  this 
reason  no  terms  involving  these  displacements  appear  inside  the 
sun  in  Eq  (A. 2).  As  a  consequence,  these  displacements  do  not 
affect  the  coherence  properties  of  the  scattered  beam  and  the 
diffraction  line  shape.  Note  that  in  its  simplest  form  this 
phenomenon  appears  in  an  interference  experiment  in  which  we 
superimpose  the  interference  patterns  generated  by  a  set  of 
screens  with  two  pinholes.  In  each  screen  the  pinholes  are  off, 
by  a  distance  fiR,  from  their  intended  position;  the  probability 
that  on  a  given  screen  the  error  6R  has  a  specified  value  is 
given  by  a  Gaussian;  there  is  no  correlation  between  the  errors 
made  for  the  pinholes  on  the  same  screen. 

The  Debye-Waller  factor  exp{-2W(^)}  appearing  in  Eq .  (A. 2) 

is  intimately  connected  to  perturbation  theory.  The  general 
phenomenon  is  not.  It  simply  has  to  do  with  the  fact  that  the 
elastic  scattering  intensity  is  depressed  because  the  the  thermal 
motion  makes  the  atomic  grattlng  slightly  faulty.  The  effect 
does  not  modify  the  line  width  because  the  errors  made  at 
different  sites  are  uncorrelated.  Thus,  we  can  generalize  the 
effect,  as  done  here,  by  using  the  Langevin  equation  to  generate 
the  random,  uncorrelated  motion  of  the  lattice  atoms  and  not 
perturbative  scattering  theory  to  describe  its  effect  on 
diffraction.  The  previous  work  followed  a  different  path.  It 
assumed  that  with  adequate  repair  the  neutron  scattering  formula 
(A. 2)  can  be  applied  to  atom  scattering.  (a)  Since  perturbation 
theory  is  not  expected  to  work  the  formula  was  corrected  in  the 
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Spirit  of  the  distorted  wave  approximation.  It  was  assumed  that 
the  atom-lattice  interaction  has  two  effects:  first  it  speeds  up 
the  incident  particle;  then  the  accelerated  particle  is  scattered 
by  the  hard  core  potential  according  to  a  formula  of  the  type 
(A. 2).  Thus  Eq.  (A. 2)  (or  its  version,  corrected  to  deal  with 
surface  rather  then  bulk  scattering)  can  be  used  if  the  initial 
momentum  appearing  in  it  takes  into  account  the  acceleration 
mentioned  -above.  This  is  chosen  to  give  a  kinetic  energy  equal 
to  the  initial  kLa-etic  energy  plus  the  well  depth.  First  it  is 
not  clear  what  permits  us  to  divide  the  collision  processes 
precisely  in  this  way.  One  may  argue  that  perhaps  we  should  take 
the  initial  momentum  zero  since  at  the  moment  of  impact  the 
particle  had  been  slowed  down  by  the  repulsive  part  of  the 
potential.  Furthermore  the  interaction  with  a  potential  of  a 
finite  range  1  does  not  create  a  planar  wave  with  a  large  wave 
vector  but  a  wave  packet  of  with  1~^  (in  momentum).  Therefore  we 
should  not  use  the  E,q .  (A. 2)  but  a  version  in  which  the  Initial 
wave  packet  is  a  Gaussian.  This  would  lead  to  different 
results.  Finally,  on  more  formal  grounds,  we  should  use  a  T 
matrix  rather  than  V  and  there  is  no  clear  reason  why  T(?)  should 
resemble  V(T)  (e.g.  should  have  a  "potential  well",  etc.)  An 
examnation  of  the  second  order  terms  <V>G<V>  in  the  expansion  of 
the  T  matrix  in  terras  of  <V>  shows  clearly  that  the  higher  order 
terras  do  not  contain  a  DW  factor  of  the  type  given  by  Eq.  (A. 2). 

b)  Since  the  pair-wise  interaction  between  the  incident 

atom  and  the  solid  atoms  is  not  as  well  localized  as  the  neutron- 

atom  one,  it  is  impossible  to  avoid  the  simultaneous  interaction 

of  the  incident  atom  with  several  solid  atoms:  moreover 

successive  interactions  with  several  atoms  are  also  important. 

None  of  these  "multiple  collision''  events  lead  to  a  Debye-Waller 

like  terra.  The  attempts  made  to  correct  for  such  effects  are  not 
2 

systematic  and  do  not  seem  to  lead  to  agreement  with  experi- 
2 

ment . 


APPENDIX  2 


In  order  to  properly  coapute  the  width  of  a  diffraction 

peak,  using  the  method  described  in  this  paper,  one  must  scatter 

a  grid  of  packets  from  a  dynamical  surface  of  sufficient  area. 

The  width  of  this  area  is  given  by  the  correlation  length  X~^, 

which  is  the  length  over  which  the  density-density  correlation 

function  for  the  surface  atoms  decays  to  the  density  squared. 

The  nature  of  this  decay  determines  the  diffraction  line  shapes. 

Making  the  primary  zone  larger  than  (X(T))  ^  is  redundant.  When 

-2  - 

we  scatter  a  bean  with  a  cross  section  smaller  than  X(T)  the 
beam  size  controls  the  diffraction  peak  width.  We  demonstrate 
this  with  the  following  simple  model  taken  from  neutron 
scattering . 

For  neutron  scattering,  first  order  perturbation  theory  is 
generally  used  to  relate  the  scattering  probability  to  the  above 

4 

mentioned  density  density  correlation  function.  One  finds. 


P(k' )  a  Jdt  e 


•lut 


1 

IJ 


<e 


+  (t) 


-lAk*Rj (0) 


where  ?.(t)  is  the  tine  dependent  position  of  the  1th  lattice 

^  ^ 

atom,  and  Ak  •  k'  -  k  is  the  change  in  wave  vector  of  the 
scattered  neutron  beam.  Writing  this  for  one  dimension,  where  ^ 
■>  nc  (n-0,l,...N)  labels  the  equilibrium  position  of  the  nth 
atom,  we  assume  that  the  above  correlation  function  for  the 
fluctuations  from  equilibrium  decays  exponentially  in  space  with 
some  length  X  Thus, 


N 

P(K)  a  Z 


^-iAk(n-m)c  ^-Xln-m)c 


n  ,  a  -  0 

where  c  is  again  the  lattice  spacing.  We  do  the  summation 
exactly  to  find 


This  Is  also  the  result  when  the  beam  size  is  smaller  than  x”  , 
(Nc<X  ^),  and  the  surface  is  coherent  as  far  as  the  probe  can 
detect.  In  this  case  the  peak  widths  are  independent  of  X{T), 
having  widths  in  k-^space  of  2(Nc)”^.  As  N  and  X~^  become  very 
large,  we  get  delta  function  peaks. 

In  the  opposite  limit,  the  beam  covers  an  area  larger  than 
the  correlation  length,  and  Nc>X~^.  In  the  neighborhood  of  a 
peak,  at  a  6k  ■=  dk  -  2Trm/c.  Eq.  (A. 3)  reduces  to 


V  2N  X 

P(q)  «  - - 5 - - 


(we  assume  that  X  ^  is  larger  than  a  few  lattice  spacings).  The 
peak  structure  is  now  independent  of  N,  and  the  width  is  X(T). 
which  increases  with  temperature.  The  shape  of  the  peak  is 
determined  by  the  manner  in  which  the  density-density  correlation 
function  depends  on  distance;  the  Lorentzlan  above  being  a 
consequence  of  the  choice  of  a  simple  exponential  decay. 

Thus,  to  generate  the  proper  temperature  dependence  in  the 
diffraction  peak  widths,  we  need  to  use  the  ADT  method  to 
propagate  a  primary  surface  larger  than  (X{T))~^,  and  probe  it 
with  a  beam  of  at  least  this  size.  In  general  this  size  can  be 
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quite  karse  at  low  temperatures,  requlrine  much  computer  tinre. 
However,  in  order  to  generate  the  proper  dependence  of  the 
elastic  peak  intensities  with  temperature  due  to  lattice  motion, 
it  is  sufficient  to  scatter  from  a  single  unit  cell.  We  can  see 
this  from  a  simple  examination  of  the  neutron  scattering  cross 
section  for  a  thermal  lattice.  This  expression  can  be  written  in 
terms  of  the  density-density  correlation  function  for  the  lattice 
atoms  mentioned  earlier.  That  is 

P  oc  /  dt  e  e^^^  ^  C(lr,t) 

-«a 

where 

GiT.t)  -  ^  Sdr  <p(‘r  -r',o)p('r  ,t)>  . 

It  is  possible  to  divide  this  expression  into  a  sum  of  elastic 
(time  independent)  and  Inelastic  (time  dependent)  terms  by  taking 

G(T,t)  =  G('r,«)  G'(^.t) 

where  G(7,«*)  is  time  independent  since 

^  2 
llm  <p ( r  ,  t  )p ( o , o ) >  -  <p> 

t-MO 

The  elastic  peak  Intensity  is  thus 

P^j  a|/d'r  e*'^^  ’ '^<p  (7)  >  I  ^  o  S  6  (Ak-x )  exp  [ -2W  (T)  ] 

•* 

T 

where  7  are  the  reciprocal  lattice  vectors  of  the  solid.  The 
Debye-Waller  factor  exp[-2W(7)l,  therefore,  is  related  to  the 
Fourier  transform  of  the  density  operator  for  a  single  lattice 
site.  All  site  to  site  time  and  spatial  correlations  lead  to 
energy  and  momentum  transfer  (l.e..  Inelastic  processes),  and  are 
included  in  G  '  (7,  t )  . 
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Fleure  Captions 


Figure  1 


Figure  2, 


Figure  3, 


Figure  4, 


Figure  5, 


Figure  6, 


Figure  7 


Velocity  distribution  of  scattered  beam  relative  to 
average  initial  velocity  (Rk^/m)  for  an  intial  grid 
of  packets  covering  NxN  cells,  where  N=1  (XXXXX), 

2  ,  3( - ),  and  4( — - - ).  Surface 

temperature  is  100  K. 

Sane  as  above  for  -N=2  case,  except  the  surface 
temperatures  are  100  K  (XXXXX)  and  400  K  (000). 

Peak  probability  vs.  AK .  AK  is  in  units  of  the 
reciprocal  lattice  vector  in  the  (1,0)  direction 
(l.e.,  AK  •  1.0  marks  the  (1,0)  peak).  The  plots 
are  for  an  initial  grid  covering  NxN  unit  cells, 

where  N  =  2,  T»=10K  ( - =^)  .  N»2,  T  =  200K  ( - ),  N  =  3, 

T=10K  (0000),  and  N»3.  T-200K  (-0-0-0-0).  Note  that 
peak  height  is  independent  of  N,  which  only  affects 
the  width. 

Peak  probabilities  P(T)  as  a  function  of  surface 
temperature,  'relative  to  static  lattice  values  P(0), 
for  the  (1,0)  (^^O)  ,  (2,0)  (000),  and  (1,1)(XXXX) 

peaks  for  He/PT(111). 

Log  of  the  relative  probabilities  of  Fig.  4  vs. 
temperature,  for  the  (1,0)  (0000)  and  (1,1)  (OO^) 
peaks . 

Peak  probability  vs.  teraperture  for  the  (0,0)  (OO^ 
(2,0)  (0000),  and  (-2,3)(XXXX)  peaks  of  Ne/Pt(lll). 

Log  of  the  probabilities  in  Fig.  6  vs.  temperature 
for  the  (0,0)  (0000)  and  (2,0)  (0  O)  peaks. 
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